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ABSTRACT 

During the last decade it has become evident that the magnetorotational instability is at the heart of the 
enhanced angular momentum transport in weakly magnetized accretion disks around neutron stars and black 
l/^ ■ holes. In this paper, we investigate the local linear stability of differentially rotating, magnetized flows and the 

f^ ' evolution of the magnetorotational instability beyond the weak-field limit. We show that, when superthermal 

^^ , toroidal fields are considered, the effects of both compressibility and magnetic tension forces, which are related 

CN ' to the curvature of toroidal field lines, should be taken fully into account. We demonstrate that the presence of 

;_( ' a strong toroidal component in the magnetic field plays a non-trivial role. When strong fields are considered, 

Qh, the strength of the toroidal magnetic field not only modifies the growth rates of the unstable modes but also de- 

termines which modes are subject to instabilities. We find that, for rotating configurations with Keplerian laws, 
the magnetorotational instability is stabilized at low wavenumbers for toroidal Alfven speeds exceeding the 
geometric mean of the sound speed and the rotational speed. For a broad range of magnetic field strengths, we 
also find that two additional distinct instabilities are present; they both appear as the result of coupling between 
the modes that become the Alfven and the slow modes in the limit of no rotation. We discuss the significance 
of our findings for the stability of cold, magnetically dominated, rotating fluids and argue that, for these sys- 
tems, the curvature of toroidal field lines cannot be neglected even when short wavelength perturbations are 
r — I considered. We also comment on the implications of our results for the vahdity of shearing box simulations in 

f^ , which superthermal toroidal fields are generated. 

^~^ • Subject headings: accretion, accretion disks — MHD — instabilities — plasmas 

o 

,-J^ ■ 1. INTRODUCTION 

Linear mode analyses provide a useful tool in gaining important insight into the relevant physical processes determining the 
^ stability of magnetized accretion flows. Studies of local linear modes of accretion disks threaded by weak magnetic fields have 
^ pointed out important clues on viable mechanisms for angul ar momentum tr anspor t and the subsequent accretion of matter 
j§ , onto flie central objects fBalhus & Hawlev" 199lL fl998l l20(32t fSano & Mivamalll999t lBalbusl l2003). They have also provided 
r • ' simplified physical m odels and analogies over which more complex physics can, in principle, be added (Balbus & Hawlev 1993, 
[1998; Ouataert. Dorland. & Hammett 2002). These treatments were carried out in the magnetohydrodynamic (MHD) limit (but 
see Quataert, Dorland & Hammett 2002 who studied the kinetic limit) and invoked a number of approximations appropriate to 
5^] the study of the evolution of short-wavelength perturbations, when weak fields are considered. In this context, the strength of the 
d magnetic field, B, is inferred by comparing the thermal pressure, P, to the magnetic pressure and is characterized by a plasma 
parameter, j3 = 8n P/B^ > 1. 

It is not hard to find situations of astrophysical interest, however, in which the condition of weak magnetic fields is not satisfied. 
A common example of such a situation is the innermost region of an accretion disk around a magnetic neutron star. It is widely 
accepted that X-ray pulsars are powered by accretion of matter onto the polar caps of magnetic neutron stars. For this to occur, 
matter in the nearly Keplerian accretion disk has to be funneled along the field lines. This suggests that, at some radius, centrifugal 
forces and thermal pressure have to be overcome by magnetic stresses, leading naturally to regions where f3 <1. 

In the context of accretion disks, the presence of superthermal fields in rarefied coronae also seems hard to avoid, if coronal 
heating is a direct consequence of the internal dynamics of the disk itself rather than being produced by external irradiation 
from the central object. Three-dimensional MHD simulations bv lMiller & Stond ( 1200 0) showed that magnetic turbulence can 
effectively couple with buoyancy to transport the magnetic energy produced by the magnetorotational instability (MRI) in weakly 
magnetized disks to create a strongly magnetized corona within a few scale heights from the disk plane. On long time scales, 
the average vertical disk structure consists of a weakly magnetized (/? ~ 50) turbulent core below two scale heights and a 
strongly magnetized {(3 < 0.1) non-turbulent corona above it. The late stages o f evolution in these models show that the disks 
themselves become magnetically dominated. 'Machida. Havashi. & Matsumotol fT^QQ^ also found that the average plasma f3 in 
disk coronae is ~ 0.1 — 1 and the volume filling factor for regions with < 0.3 is up to 0.1. Even in the absence of an initial 
toroidal component, simulations carried out by Kudoh, Matsumoto, & Shibata ( 2002) showed that low-/3 regions develop near the 
equator of the disk because of a strong toroidal component of the magnetic field generated by shear On more theoretical grounds, 
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strong toroidal magnetic fields produced by strong shear in the bo undary layer region have been suggested as responsible for the 
observed bipolar outflows in young stellar objects (PringlqQ^l^)- More recently, Pariev. Blackman. & Boldvrev (2003) found 
self-consistent solutions for thin magnetically-supported accretion disks and pointed out the necessity of assessing the stability 
properties of such configurations. 

Another case in which magnetic fields seem to play an important dynamical role in rotating fluid configurations is that of 
magnetically supported molecular clouds. Observations of both large Zeeman line-splitting and of broad molecular lines support 
the presence of superthermal fields (see Myers & Goodman 1988 for further references and Bourke & Goodman 2003 for a review 
on the current understanding on the role of magnetic fields in molecular clouds). Values of the plas ma of the order of 0. 1 — 0.0 1 
have also been used in numerical studies of the structural properties of giant molecular clouds (.Ostriker. Stone. & Gammiel2001h . 

As a last example of astrophysical interest, we mention magnetocentrifugally driven winds, such as those observed in protostars. 
These outflows seem to play an important role in the evolution of young stellar objects and in the dynamics of the parent clouds 
by providing a source of turbulent energy. Magnetocentrifugal jets typically involve internal Alfven speeds comparable to the 
flow speeds. These structures are supported mainly by magnetic pressure due to strong toroidal fields. The ratio of magnetic 
pressure in the jet to the gas pressure of the ambient mediu m ca n be of the order of 10^ (for an extensive study of MHD driven 
instabilities in these systems see Kim & Ostriker 2000 and ii5.3> . 



In this paper, we investigate the local linear stability of differentially rotating flows without imposing any a priori restrictions 
on the strength of the magnetic field. We do, however, restrict our attention to rotationally supported flows (we loosely use this 
term to refer to flows with internal Alfven speeds smaller than the rotational speed). Our intent is to demonstrate that the effects of 
the finite curvature of the toroidal field lines on the stability of small-wavelength vertical perturbations (i.e., on the most unstable 
modes present in the weak-field MRI) cannot be neglected when superthermal toroidal fields are present. In order to achieve 
this task, we relax the Boussinesq approximation (see also Papaloizo u & Szuszkiewicz 1992 and Blaes & Balbus 1994), which 
is valid only when the toroidal component of the field is subthermal (Balbus & Hawlev 1991). We thus consider the MHD fluid 
to be fully compressible. Moreover, even though we perform a local analysis, we do consider curvature terms when evaluating 
magnetic forces, for they become important in the strong-field regime (see also Knobloch 1992 and Kim & Ostriker 2000). 

In most early studies addressing the MRI, it was found that the only role played by a toroidal component in the magnetic field 
is to quench the growth rates of the modes that are already unstable when only weak vertical fields are considered (Balbus & 
Hawley 1991; Blaes & Balbus 1994; see Quataert, Dorland, & Hammett 2002 for the kinetic limit; see also Kim & Ostriker 2000 
for the cold MHD limit). Here, we show that, when strong fields are considered and the approximations usually invoked in the 
study of the weak-field MRI are relaxed, the presence of a toroidal component of the magnetic field plays a crucial role not only 
in the growth rates of the unstable modes but also in determining which modes are subject to instabilities'*. As expected, the 
presence of a toroidal component breaks the symmetry of the problem, also giving rise to traveling modes. Moreover, for a broad 
range of magnetic field strengths, we find that two different instabilities are present. They both appear as the result of coupling 
between the modes that become the Alfven and the slow mode in the limit of no rotation. 

The paper is organized as follows. In fJ2l we describe the physical setup to be studied, present the dispersion relation to 
be solved, and discuss the importance of curvature terms in the limit of superthermal fields. In fJ3] we solve numerically the 
dispersion relation in some interesting regimes. In fJ4] we study the onset of instabilities as a function of magnetic field strengths 
and present some useful approximate criteria that enable us to study analytically some aspects of the full problem. In fJ5] we 
compare our results to previous investigations and discuss some of the implications of this study. Finally, in %\ we present a brief 
summary and our conclusions. 

2. MHD EQUATIONS FOR PERTURBATIONS AND THE DISPERSION RELATION 



We Start with the set of equations that govern the behavior of a polytropic fluid in the MHD approximation 

dp 

dt 



^ + V-{pv) = 0, (1) 



dv , _, _, _/„ B'^\ (B 



P- + (P..V) . ^ -pV$ - V ^P + ^ j + ^^-V j B , (2) 

^- + (V-u) B - (B- V) -u + (-U-V) B = , (3) 

at 

and 

^-'"{il ■ <^> 

In these equations, p is the mass density, v the velocity, P the gas pressure, and F the polytropic index; B is the magnetic field 
and $ the gravitational potential. For convenience we adopt a cylindrical set of coordinates (r, 0, z) with origin in the central 
object (i.e., the neutron star or black hole). We assume a steady axisymmetric background flow characterized by a velocity 

^ We will comment later in more detail on the paper bv lCurrv & Pudrita I1995l> who outlined the effects of a dynamically important toroidal field (in the case 
of an incompressible MHD flow) and address the similitudes and differences with our findings. 
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field of the form v = v^{r, z)(/) and threaded by a background magnetic field. For consistency, our analysis is restricted to 
background fields of the form B = Ba<J) + B^z since the effect of including a radial component in the field is to generate 



a linear growth in time of the toroidal component (Balbus & Hawlev 1991). Under these circumstances, all the background 
quantities depend on the radial and vertical coordinates only. In the present treatment, we neglect the self gravity of the fluid. In 
fact, |Pariev, Blackman, & Boldvrev (2003) showed that magnetically dominated accretion disks have lower surface and volume 
densities for a fixed accretion rate. This suggests that these systems are lighter than standard disks and thus are not subject to 
self-gravity instabilities. 



2.1. Equations for the Perturbations 

In order to perform the local linear mode analysis, we perturb the set of equations Q-Q by substituting every physical variable 
fhyf + Sf and retain only linear orders in Sf. In the following, we focus our analysis on the study of axisymmetric perturbations 
in an axisymmetric background. 

We first consider, in some detail, the radial component of the momentum equation (|2}, which becomes 
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(5) 



The coefficients in this linear equation for the perturbed variables depend in general on r and z, [e.g., the angular velocity Vl 
is in general r2(r, z)]. Therefore, at this point, the decomposition of perturbed quantities in Fourier modes - e.g., symbolically 
(5/ = ^ ^j'^g«(feri'+/c^z-tjt) _ y^ould not result in any particular simplification of the problem. This can be seen by taking the 
Fourier transform of equation (|5} which results in a sum of convolutions of the Fourier transforms of background and perturbed 
quantities. Further progress can be made if we restrict the wavelengths of the perturbations for which our stability analysis is 
valid. To this end, we choose a fiducial point, ro — (ro, (t>o, zq), around which we perform the local stability analysis. The choice 
of the particular value of 00 is, of course, irrelevant in the axisymmetric case under study. 

We expand all the background quantities in equation (|5} in Taylor series around ro and retain only the zeroth order in terms of 



the local coordinates S,r 

dSVr 



ro and £_, 



dt 



2VIq5va 



1 

Airpo 



B' 



1 

Po 
dB, 



dSP 



dr 



B' 



1 

47r \ dr 

dB 



zq to obtain 
dB^ 



5Ba + B 



, dSB^ dB, 



dr 



5p 
Pa 



d£,r dr 

odSBr 



SB, + B\ 



, dSB, 



dP 
dr 



Anpo 



B 



d^. 



- 2^SB^ - ^ 
ro I inpo 



1 (Birsp 



ro Po 



S_p 

QPO 



= 0. 



(6) 



Here, ilo, po, B'^, and B'^ stand for the angular velocity, background density, and magnetic field components at the fiducial point 
ro and the subscript "0" in the derivatives with respect to the radial coordinate r indicates that they are evaluated at ro. Equation 
(|6jl is a linear partial-differential equation in the local variables for the perturbed quantities but with constant coefficients. This is 
a good approximation as long as the departures (^,. , ^^ ) are small compared to the length scales over which there are significant 
variations in the background quantities, i.e, £_r ^ Lr and (,, <C L,, where Lr and L, are the characteristic length scales in the 
radial and vertical directions, respectively. 

It is only now that it is useful to expand the perturbed quantities in equation (|6} in Fourier modes. We can thus write each one 
of the perturbed quantities as 



■r J "^2; , t 



Sf = Sf{h, 
and write the radial momentum equation for each mode as 
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Here, we have introduced the quantities Cg , v\j, and uj(^ that stand for the local sound speed and local Alfven speeds associated 
with the toroidal and vertical components of the local magnetic field and are defined by 
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(9) 



Note that, for brevity, we have omitted the dependences (fcr, kz,u!) in the Fourier amplitudes. For consistency, the validity of 
the analysis is now restricted to modes with wavenumbers satisfying fc^ir ^ 1 and k^L, ^ 1. Without loss of generality. 
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we assume that the fiducial point is inside the disk so we can write the local conditions on the wavenumbers as krrQ 3> 1 and 
kzZQ :^ 1 . Moreover, for fiducial points such that ro > zq the latter condition also implies fc^ro ^ 1. 

At this point, it is also convenient to define a new set of independent variables (JwAr, <5wa0, <^waz) defined in terms of 
{SBr, 5B^, SBz) in such a way that 6va = SB/^Anpo. In this case, equation ^ reads 
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As a last step, it is useful to work with dimensionless quantities. To this end, we define dimensionless variables by scaling all 
the frequencies with the local rotational frequency Hq and all speeds with the local circular velocity fio'^o- It is also convenient to 
define dimensionless wavenumbers by multiplying the physical wavenumberby the radial coordinate tq. In summary, we define 



ti) ~ lo/Uq , 



fhrp Kf^lQ . 
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dp = Sp/po , Sv = Sv/floro , 6va = SvA/^oro . 

For completeness, we introduce here the epicyclic frequency k and its local dimensionless counterpart 
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(12) 
(13) 



(14) 



This quantity appears naturally in stability analyses of differentially rotating configurations and it is the frequency at which all 
the flow variables oscillate around their background values in the absence of magnetic fields. For a rotational profile given by a 
power law, the epicyclic frequency is proportional to the angular frequency at all radii. 

Finally, the dimensionless version of equation ( I10> reads. 
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Following a similar procedure with the remaining equations in the system Q-(0}, we arrive to the linear set of equations needed 
to perform the local stability analysis. For brevity, we now drop the hat in all the dimensionless variables and the superscript in 
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where we have used equation to recast the pressure perturbations in terms of density perturbations. The local conditions over 
the wavenumbers now read kr ^ 1 and k^ ^ I. 

The factors e,, with i — 1,2,3,4, are just convenient dummy variables that we introduce in order to help us keep track of the 
terms that account for the finite curvature of the background and are usually neglected in local studies of the weak-field MRI. 
Their numerical values are to be regarded as unity, unless otherwise mentioned. The terms proportional to ei and 62 in equation 
( I17> and the term proportional to £3 in equation (I18> are due to the effects of magnetic tension and they appear naturally when a 
cylindrical coordinate system is adopted. The terms proportional to £4 in equations (I16> and (I22> are related to flux conservation 
in cylindrical coordinates. Although the three terms labeled by ei, e2, and £3 share the same physical origin (i.e., magnetic tension 
introduced by the curvature of toroidal field lines), it is useful to be able to distinguish among them because the one labeled with 
£2 vanishes in the limit of an incompressible flow. Note that equations (|20} and (12 2 > ensure a divergence-free perturbed magnetic 
field, i.e., 'V-SB — 0, only when the finite curvature of the background is accounted for (i.e., £4 = 1). 

Up to this point, our intention has been to keep the discussion as general as possible in order to clearly state all the assumptions 
that we have made to obtain the set of equations for the perturbations to perform a local linear mode analysis. For the sake of 
simplicity, and to avoid the parametric study from being too extensive, we further invoke the following assumptions. We choose 
the fiducial point Tq to lie in the disk mid-plane and assume that, locally, the vertical gradients in all background quantities are 
negligible and set them to zero. This will be a good approximation as long as we consider equilibrium configurations such that 
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Note that the solenoidal character of the magnetic field ensures that, for B = B^{r, z)(f) + Bz{r, z)z, the condition dB^/dz — 
holds for arbitrary z. 

In general, the forces induced by the curvature terms (e.g., the one proportional to £2 in eq. 1171 ) and those induced by 
background (logarithmic) gradients in the radial direction, (e.g., the term proportional to dlnS^/dlnrjo in the same equation) 
will not cancel each other. As a single exception, for the case in which B^p oc r^^, the most important effects due to the finite 
curvature of toroidal field lines are canceled out by the gradients in the toroidal field. This particular case, however, might not be 
completely relevant to realistic rotating flows since, in order to ensure force balance when the thermal pressure can be neglected 
against magnetic stresses, the magnetic field strength must decline more slowly than r^^ (see, e.g., Kim & Ostriker 2000). For 
simplicity, we further focus our attention on the study of differentially rotating, axisymmetric MHD flows with locally negligible 
radial gradients in the background density and magnetic field, i.e.. 
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In the rest of the paper, we consider that the only background flow variable with a non-negligible local radial gradient is the 
angular velocity Vt ex r^'^ and set all other radial gradients to zero. Note that this assumption is widely invoked in many 
investigations of the weak-field MRI (e.g., Blaes & Balbus 1994; Balbus & Hawley 1998; Blaes & Socrates 2001; Quataert, 
Dorland, & Hammett 2002; Balbus 2003). The assumption that the only background variable with a non-negligible radial 
gradient is the angular velocity is also generally a part of the initial set of conditions used in many numerical analyses of the MRI 
in the shearing box approximation (e.g., Hawley, Gammie, & Balbus 1994, 1995, 1996; Miller & Stone 2000). 

In spite of being linear in the perturbed quantities, the terms proportional to €i have been neglected in previous local studies 
of the MRI under the assumption that fc^ ^ 1 and k^ ':> 1 (but see also Knobloch 1992 and Kim & Ostriker 2000). Although 
comparing an imaginary term against a real one in a stability analysis might seem particularly risky, this might not be a bad 
argument in order to neglect the terms proportional to £1 in equation ( I17> or £4 in equations ( I16> and (12 2> against ikr (but see 
the discussion in Appendix A). The same could be said about the terms proportional to £2 in equation ( I17> or £3 in equation 
( I18> in the limit of a very weak toroidal component in the magnetic field, given that both of them are proportional to va4,- It is 
not evident, however, that we can neglect the terms proportional to either £2 in equation ( I17> or £3 in equation ( I18> if we are to 
explore the regime of strong toroidal fields. There are two different reasons for this. In order to neglect the term proportional 
to £2 against the one proportional to fc^ in equation (^} we should be able to ensure that the condition {^2lkr){v\./c^) ^ 1 is 
always satisfied, since both terms are proportional to 6 p. In this particular case, neglecting the forces induced by the bending of 
toroidal field lines becomes a progressively worse approximation the colder the disk is and is not well justified in the limit Cg -^ 0. 
The case presented in equation (I18> is even harder to justify a priori since now we would need to guarantee that the condition 
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{£3/kz)ivAct,/vAz){SvAr/SvA<t>) ^ 1 is always satisfied. However, this ratio is not only proportional to va4,/vaz, which might 
not be negligible in many astrophysical contexts but, through the ratio Svat/Svac/i, is also a function of kr, kz, and uj{kr, kz); the 
magnitude of this term is therefore unknown until we solve the problem fully. A similar situation to th is one is encountered if we 
aim to compare the term e2v\,5p with the term proportional to kzVAzSvAr in equation ( I17> (see ilS.ll for further discussion). 

For the sake of consistency and in order not to impose a constraint on the magnitude of the toroidal Alfven speed with respect 
to the sound speed we keep all the terms proportional to the parameters e^. We will later show that the term proportional to ei is 
negligible when superthermal toroidal fields are considered. We will also discuss under which conditions the terms proportional 
to £4 can be neglected and why the terms proportional to €2 and e^ are particularly important. 

2.2. Dispersion Relation 

In order to seek for non-trivial solutions of the homogeneous system of linear equations ( I16> - (I22> we set its determinant to 
zero. The resulting characteristic polynomial is 
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where we have dropped the subscript "0" in the radial logarithmic derivative of the angular frequency. This is the most general 
dispersion relation under our current set of assumptions. When all the parameters Ci are set equal to zero, we recover the results 
of previous analyses where the curvature of the toroidal field lines was not considered (e.g., Blaes & Balbus 1994; Balbus & 
Hawley 1998), while when they are set equal to unity we obtain our full dispersion relation. 

Although the original hnear system ( I16t - (l22t related seven variables (recall that we had eliminated 5P in terms of 5p using 
eq. which is time-independent), the characteristic polynomial is only of 6'^ degree. This is easily understood by noting that 
equations MQ\ and ( I22> can be combined into one single equation expressing the solenoidal character of the perturbations in the 
magnetic field, 'V-6B = 0. This implies a relationship between SBr and SBz (or equivalently between SvAr and 5vaz) that 
must be satisfied at all times and is, therefore, independent of uj. The fact that the dispersion relation ( I25l l is of 6*'' and not of 
4'^ degree is because we are taking into account the effects of finite compressibility. This can be seen immediately by taking the 
limit Cs ^ 00. 

Once all the dimensionless variables have been properly defined, it is not evident that the magnetic -tension terms, proportional 
to ei, 62, and 63, will play a negligible role in determining the eigenfrequencies uj. This is because the non-vanishing toroidal 
component of the magnetic field introduces odd powers in the dispersion relation and hence break its even symmetry. In fact, 
small modifications in the odd-power coefficients can and do have an important impact on the nature (real vs. complex) of the 
solutions. As we will see in fjsjand describe in further detail in fj3 these curvature terms introduce further coupling between the 
radial and toroidal directions, which in turn result in a strong coupling between the Alfven and the slow mode. 

Also important is the fact that some of the coefficients in the dispersion relation ( I25> are no longer real due to the factors ikr- 
The presence of these terms does not allow us to affirm that complex roots will appear in conjugate pairs. As we discuss in 
Appendix A, the terms proportional to ik,- play an important role in determining the stability of modes for which the ratio kr/kz 
is non-negligible, even in the local limit, i.e., when kr ^ 1. Of course the smaller the ratio kr/kz, the smaller the effects of the 
factors ikr will be. If we consider the limit kz ^ kr in equation ( I25> . the imaginary part of all the coefficients in the dispersion 
relation will become negligible. In this limiting case, whenever a given complex root is a solution of the dispersion relation ( I25> 
so is its complex conjugate, for the dispersion relation has real coefficients (see Appendix A). 



In the next section, we will show that the dispersion relation ( I25t reduces to the dispersion relations previously derived in many 
local studies in different regimes. It is important to emphasize, however, that this dispersion relation fully considers the effects 
of compressibility and magnetic tension simultaneously without imposing any restrictions on the field strength or geometry. This 
feature is crucial in determining the stability properties of the MHD flow when strong toroidal fields are considered. 

2.3. Previous Treatments 

There has been some discussion in the past abo ut the importance of the curvature terms for the stability of magnetized Ke- 
plerian flow s (Knobloch 1992; Gammie & Balbusl 1 1994ft . In studies in which these terms were considered (Knobloch 1993 
iDubrufle & Knobloch.199 3). compressibility effects were neglected. On the other hand, t here h ave also been treatments in which 
compressibility was addressed but the curvature terms were neglected (Blaes & Balbus: [l994h . Both types of studies provided 
arguments for and against the importance of these terms. The limit of cold MHD flows has been addressed by Kim & Ostriker 
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(2000). These authors concluded that when the magnetic field strength is superthermal, the inclusion of toroidal fields tends to 
suppress the growth of the MRI and that for quasi-toroidal field configurations no axisymmetric MRI takes place in the limit 

Because of the generality of our treatment, in which both curvature terms and compressibility effects are fully taken into 
account, we are able to address all of these issues in SjS] For the time being, and as a check, we can take the appropriate limits in 
the general dispersion relation ( I25> to recover the dispersion relations derived in the aforementioned works. 

Compressibility with no field curvature — Setting e^ — 0, for i = 1,2, 3, 4, and considering perturbations propagating only in 
the vertical direction (this can be formally done by taking the limit fc^ ^ kr in equation 1251 ) we recover the dispersion relation 
derived in the compressible, weak-field limit bv lBlaes & BalbusI J 19941) . 

co'-[kl{c'.+vi, + 2vlJ + n']u' 



d\n.r 



u? 



,dlnf7 



fe'^LcsM fc>L + 2^^^ I = . (26) 



The stability criterion derived from this dispersion relation is not different from the one derived, within the Boussinesq approxi- 
mation, by Balbus & Haw lev ( 1991 ). All the perturbations with vertical wavenumber smaller than the critical wavenumber /cbh 
are unstable, with 






kln^L ^ -2^^^ . (27) 



In this case, the strength of the toroidal component of the magnetic field does not play any role in deciding which modes are 
subject to instabilities. 

Field curvature with no compressibility — It is important to stress that even in the incompressible limit not all the terms 
proportional to e^ in the dispersion relation ( I25> are negligible (of course, the ones proportional to €2 are). To see that this is the 
case, we can take the limit Cg ^ 00 in the dispersion relation ( I25> to obtain 

{k1 + k'^ ~ ikreiJLU* — kl[2v'X^{kl + k'^ — e^ikr) + k^ + 2eieiv\^]uj'^ — 2fcfwA0'yAz(2ei + £3 + <^a)^ 



+k^vp^^ 



{kl + kl- eiikr)v\^_ + 2-^j^ - 261631;^^ 



, (28) 



where we have explicitly left the factors e^ that should be considered as unity. This incompressible version of our dispersion 
relation is to be compared with the one obtained by Dubrulle & Knobloch (1 993 ) as the local limit of the co rresponding eigenvalue 
proble m. Note that, in order to compare expression ( I28t with the dispersion relation (eq. [37]) presented inlDubrulle & Knoblochl 
i 19931) . it is necessary to considered the limit dv^^/dr, dvp^^/dr — > in their equation (9). We also note that the radial 
wavenumber kj. appears in equation ( I28l l only in the combination kl + k^ — ikj. while in equation (37) in fDubruUe & Knoblochl 
tl993.) we only find it as k^ + k"^ (i.e., ?i^ + k"^ in their notation). This is because when taking the local limit, kr 3> 1, in 



the process of deriving their equation (37) from their equation (9), the terms proportional to ikj. were neglected against fc^ by 
iDuhmlle&KnobloclJJI 993ft. 



When the toroidal component of the magnetic field is negligible, i.e., when vx^j, -^ in equation ( I28> . and we consider vertical 
modes {k^ ^ kr), we recover the dispersion relation for the incompressible MRI; the onset of unstable modes is still given by 
expression Mil . For weak toroidal fields, i.e., when VA<f, ^ 1, we can read off the small corrections to the critical wavenumber 
from the constant coefficient, 

{kffvl^ = -2^ + 2e,e,v\, . (29) 

For stronger fields, however, the oj — Q mode is no longer unstable (see the Appendix B for a general discussion on the stability 
of the a; = mode when compressibility and curvature terms are considered) and it is necessary to solve equation ( I28> in order to 
find the critical wavenumber for the onset of the instability. Roughly speaking, we would expect the solutions of equation (I28> to 
depart significantly from the solutions to the incompressible version of the dispersion relation \26\ when v\, > |(ilnil/(ilnr|, 
or WA0 Jc^ 1-2 for a Keplerian disk. Since in this paper we consider rotationally supported configurations (i.e., VA4, ^ 1), we will 
not address the modifications to the mode structure caused by curvature terms in incompressible MHD flows. 

It is important to stress that, for both dispersion relations ( I26> and ( I28> . in the case of rotationally supported disks, the stability 
criterion is insensitive (or, at most, very weakly sensitive, in the incompressible case) to the magnitude of the toroidal component 
of the field. As we will see throughout our study, the stability criteria that emerge from equation ( I25> are significantly different 
from the ones discussed in this section, when we consider fields for which VA4, > Cg. We will also see that the term proportional 
to €2, which depends on curvature and compressibility effects and is, therefore, absent from either equation J26> or (I28> . plays an 
important role in determining the mode structure in the general case. 

Cold limit with no field curvature — Another limit of interest is the one corresponding to the cold, MHD, cylindrical shearing 
flows usually involved in the modeling of cold disk winds (i.e., far away from the disk). In this context, Kim & Ostriker (2000) 



8 



PESSAH & PSALTIS 



addressed the behavior of the compressible axisymmetric MRI in the limit Cg -^ 0. These authors obtained a dispersion relation 
considering both vertical and radial wavenumbers and derived the criterion for instability associated with it. Their dispersion 
relation in the fully compressible case [eq. (57)] reads 



UJ 



^■']lo^ 



{kl + kDvlA'icl + v% + 4J + K^icl + vl^) + 2 
dhiQ. 



d\un 

rflnr 



VAz 



iK + KX. + 2 



d\nr 



= 0, 



(30) 



This dispersion relation can be obtained from equation (I25> if we set e^ ~ 0, for i — 1,2,3,4. Note that if we take the limit 
kz ^ kr in equation ( I30> we recover equation j26> . 

For extremely cold flows we can take the limit Cg ^ in equation ( I30> to obtain 



W* - [{kl + kl){vl^ + vl,) + fc24^ + k2]^2 ^ fc2^|jfc2 ^ ;,2)(^|^ _^ ^IJ ^ K^fc^^ii^ + 2fc2i,| 



d\nn 



dlnr 



= 



(31) 



On the other hand, taking the limit Cg 

■6 {(fc2 + fc,2)(^2 



W 



in equation (I25I I. we obtain the more general dispersion relation 

■ «L) - ifcr[(2ei - e2)v\^ + e4wL] + ^^^L + ^^ + e2e4^'i0}w^ 
(2ei + £3) 2k^VA4,VAzUj''^ + S fc^^'lzil^f + fc^ - eiikr){v\^ + wi^) + ifcr(e2 - ea)^'!^] 



2 2 



'2/cfwA0fAz(e2 



e2e4WA 



dlnr 

„4„,2 



«Az + (£264 - 2eie3)wA0f A2 



e3)wA0'^ + k^Vj^^e2e3Vj^^ = 



(32) 



Note that in this expression, as it was also the case in the incompressible limit, several of the terms that are due to the finite 
curvature of the toroidal field lines are still present. 

Analyzing the limit Cg ^ in the dispersion relation ( I30> . Kim & Ostriker (2000) concluded that toroidal fields tend to suppress 
the growth of the MRI and that, for a Keplerian rotation law, no axisymmetric MRI occurs if i < 30°, where i is the local pitch 
angle of the magnetic fields defined by i = tan~^(uAz/wA0)- However, the eigenfrequencies satisfying the dispersion relations 
j25> and ( I30> in the limit Cg — > are different and so are the criteria for instability which they are subject to. In MS. 31 we comment 
in more detail on how the solutions to the dispersion relations ( I25> and ( I30> differ in the limit Cg ^ and on the implications 
regarding the stabilization of the MRI in cold MHD shearing flows. 

In order to investigate how previous results from local stability analyses of the weak field MRI are modified as the strength 
of the toroidal field component increases, we will focus our attention on the stability of modes with k^ ^ kr^ This approach 
is physically motivated, since vertical modes correspond to the most unstable modes in the well studied MRI, and is also more 
tractable mathematically. In the next two sections, we will perform a thorough numerical and semi-analytical study of the general 
dispersion relation ( I25> in the limit kz ^ kr, with particular emphasis on the case of strong toroidal fields. We will then be in a 
better position to understand the similitudes and differences of our findings with those of the aforementioned studies and we will 
address them in |5] 



3. NUMERICAL SOLUTIONS 

We solved numerically the dispersion relation ( I25> for the frequency uj as a function of the wavenumber k^, employing La- 
guerre's root finding method (Press et al. 1992). As a typical situation of interest, we consider a Keplerian disk with Cg = 0.05 and 
VAz = 0.01. As it will be seen from the range of values of kz in which the various instabilities occur, the case of quasi-toroidal 
superthermal fields is perfectly suited to be studied in the local approximation, i.e., when k^ ^ 1, provided that the vertical 
component of the magnetic field is weak enough (i.e., vaz <C 1). 

To better appreciate the effects that the curvature terms have on the stability of the modes, a set of solutions to the dispersion 
relation j25> is shown in Figures [2 and |2^ Each of the three panels, in both figures, shows the real and imaginary parts of the 
solutions for different values of the toroidal field strength, parameterized by ua0- The left panel shows the solutions to the full 
dispersion relation ( I25> . i.e., when ei = £2 = £3 = £4 = 1- The central panel shows the solutions to equation ( I25> when 
compressibility is neglected in the curvature terms, i.e., when ci = e^ — e^ = 1 and £2 = 0. For the sake of comparison, the 
right panel shows the solutions to the dispersion relation ( I26> . in which all curvature terms are neglected. 

We first analyze Figure[2 When all magnetic tension terms are neglected (right panel), the qualitative structure of the normal 
modes of the plasma is insensitive to the magnitude of the toroidal field component (see Blaes & Balbus 1994). However, the 



* In Appendix A, we briefly describe how this results are modified when finite ratios kr/kz are considered. 

' Some animations of the results presented inPigs. lTlllllTl andlslare available at http : //www. physics . arizona .edu/'- 
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Fig. 1. — The real parts of the numerical solutions to the dispersion relation 1251 coiTesponding to a Keplerian disk with Cg = 0.05 and tiAz = 0.01. Leftpanel: 
solutions to the full problem (ei = £2 = ^3 = ^4 = !)■ Central panel: the case in which compressibility is neglected in the curvature terms (ei = €3 = £4 = 1 
and £2 = 0). Right panel: the case in which all curvature terms are neglected {e\ = €2 = ^3 = ^4 = 0). Open circles indicate unstable modes (i.e., those with 
positive imaginary part). Long-dashed, short-dashed, and point-dashed lines show the fast, Alfven, and slow modes, respectively, in the limit of no rotation. 



situation is very different when the magnetic tension terms are included. For weak toroidal fields, i.e., when wa^ ^ 0.1, the 
solutions seem quite insensitive to the curvature terms; indeed these terms do not seem to play a significant role in altering the 
local stability properties of magnetized Keplerian flows compared to what is quoted elsewhere in the literature. As we will see 
later, for a Keplerian disk, the presence of the curvature terms is significant once v\j^ > Cg, which in this case translates into 
VA^ > 0.22. 
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Fig. 2. — The imaginary parts for the cases discussed in FigurelTI 



For stronger toroidal fields, i.e., when VA<f> > 0.2, the modes with the longest wavelengths become stable when all curvature 
terms are included, in sharp contrast to the case in which €2 = 0. For even stronger toroidal fields, i.e., when wa0 ^ 0.3, a second 
instability appears at long wavelengths, while the original instability is suppressed. When VAcf, J^ 0.4, both instabilities coexist as 
separate entities and the original instability reaches smaller and smaller spatial scales, when the magnitude of the toroidal field 
increases. For even higher toroidal fields, i.e., when v^^ > 0.7, the largest unstable wavenumber of the instability that developed 
for VA4, ^ 0.3 approaches fcen (see eq. 1271 ). The major implication of neglecting compressibility in the curvature terms is that 
the original instability seems to be totally suppressed for toroidal fields larger than the ones corresponding to VAcft ^ 0.3. 
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As it is clear from the dispersion relation ( I25t . the presence of the toroidal component in the field introduces odd powers of 
the mode frequency lo and hence breaks the symmetry between positive and negative real parts of the solutions. The physical 
meaning of this is clear The phase velocities of the instabilities are no longer zero and they are propagating vertically throughout 
the disk. This, of course, is not the case for the unstable solutions to the dispersion relation (I26> regardless of the magnitude of 
UA0- In that case, the most noticeable effect of an increasing toroidal field is to reduce the phase velocity of the stable modes 
beyond /cbh (which is itself independent of ua^)- 

It is also interesting to analyze how the presence of the curvature terms modifies the growth rates of the unstable modes as 
a function of the toroidal magnetic field. This is shown in Figure |2] Again, there are no significant changes for ua0 ;$ 0.1; 
however, quite significant modifications to the growth rates are present for va^ }t 0.2. The sequence of plots in the left panel 
shows more clearly the suppression of the original instability, the appearance of the instability at low wavenumbers, the return 
of the instability at high wavenumbers, and finally the fusion of these last two. The right panel in this figure shows the effects 
that the presence of a strong toroidal component has on the mode structure when curvature terms are not considered. In this 
case, the critical wavenumber for the onset of instabilities is not modified while there is a clear reduction in the growth rate of 
the non-propagating unstable modes as the magnitude of the toroidal field component increases. When the curvature terms are 
considered fully, the effects are more dramatic. Note also that the growth rate of the original instability is reduced faster from the 
first to the second plot in the left panel in Figure|2lwith respect to their counterparts in the right panel of the same figure. 

4. THE ONSET OF INSTABILITIES 

4.1. Unstable Modes 

In fJ3lwe presented how the structure of the various modes evolves as a function of the toroidal field strength and noted that, for 
a range of field strengths, two different instabilities are clearly distinguishable. Here, we obtain the conditions (i.e., the range of 
wavenumbers and toroidal field strengths) for which these unstable modes are present. We start by plotting in Figure|3]the range 
of unstable wavenumbers as a function of the toroidal field strength. As a reference, we have plotted the case for a Keplerian 
disk. The black dots in the diagram represent the unstable vertical wavenumbers, in units of k^Vj^z/cs, for a given toroidal Alfven 
speed, in terms of va<p/cs- Three regions of unstable modes are clearly distinguishable: 

• Region I shows the evolution of the original instability present in the topmost three plots in the left panel in FigurefO This 
is the region where the MRI lives. Strictly speaking, the MRI is confined to the region where ua^/cs ^ 1- As we will 
comment in 8 ^4.21 instability I is no longer incompressible beyond this point. The maximum wavenumber for which this 
instability exists is independent of VA4, and corresponds to the critical wavenumber for the onset of the MRI (i.e., /cbh in 
eq. 1271 ). The stabilization of the long-wavelength perturbations beyond a critical value of the toroidal Alfven speed is 
also evident in this region. For larger toroidal field strengths, shorter and shorter wavelengths are stabihzed up to the ones 
corresponding to fceH- 

• Region II represents the evolution of the instability that is only present for wavenumbers k^ > fceH- Note that /cbh is now 
the minimum wavenumber for the onset of instability II. In this case, increasing wa^/cs gives rise to unstable modes with 
even shorter wavelengths (two bottommost plots in the left panel of either Figure^orl^J. 

• Region III shows the instability that appears for intermediate wavenumbers (see for example the third plot in the left panel 
in Figure [0. Note that the shortest unstable wavelength in this region approaches fceH for large values of wa^/cs (i.e., 
bottommost plot in the left panel of either FigureHorl^J. 

4.2. Analytic Approximations 

In this section we obtain analytical approximations to the dispersion relation ( I25> in various limits, which will help us identify 
the different critical curves in Figure |3] 

The fast (or magnetosonic) modes are reasonably well decoupled from the rest of the oscillations (see left panels in FigQ. 
By studying the modes that satisfy the condition w^ ^ ^2^3, we effectively eliminate the fast modes from our analysis. This 
can be done for strong toroidal fields because, even in the presence of rotation, the magnetosonic modes are well described by 
cj^ ~ fcz(cs + ''^a)- Note that imposing uj'^ ^ fc^Cg is a distinct and weaker condition than asking for the MHD fluid to be 
incompressible (cg -^ 00). By eliminating these fastest modes, it is possible to find a 4*'^ degree dispersion relation in uj, with 
solutions that constitute a very good approximation to the interesting modes seen in Figures^and|2l 

We first write the equations for the evolution of the perturbations in the magnetic field. For the sake of clarity, we present the 
intermediate steps with the appropriate physical dimensions but we drop the index indicating local values. Substituting equations 
( I16> . ( I20> . and J22> in equation (I19> we obtain Sv^ in terms of 5B^ and 6Bz, 

^'~ {Kc^f-u:^ cl B^ + Bz ■ ^''^ 
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Fig. 3. — The black dots represent unstable modes obtained from solving the dispersion relation 125 i numerically for a Keplerian disk with Cs = 0.05 and 
'VAz = 0.01. The solid (fc^) and dashed (fe^*) lines con'espond to the critical wavenumber for which the a; = mode exists in the case of a compressible 
(see Appendix B) and an incompressible (discussed in i|2j flow, respectively. For strong toroidal fields, compressibility plays a crucial role in the stability of the 
LO = mode. Note that, in the limit of small vx^/cs we have fe^, fe^' — > ^bh, and the trivial mode becomes unstable. 



Using this result in equation (I21> we find. 



ik^BJv^ = - 



dlnr 



SBr — iloSBa, — iujBa- 



ykzCs) 



[k^c^Y 



■ LO^ 



From equations ( I16> . ( I22> . and ( I33> we can recast 8p in terms of 6B^ and SB^ as 



bp_ 
P 






SB, 



ihcsy 



w2 B. 



1- 



("^zCs 



(fS^Cg 



^A0 SB^ 



5. 



SB, 
B, 



Finally, we can write equations ( I33> - (I35> for the modes with frequencies such that oj'^ ^ k^c^ as. 



Sv, 



to 
kl 



(^)' 



ik^BJvA 



and 



dfl 64 uj B^ 
dlnr r kz Bz 

S_p 
P 



SB^ , . 


£4 SBr 

kzT Bz 


7 


SBr - 


l+f^ 


f). 



iloSBa 



'^Aj. SB^ 



(34) 



(35) 



(36) 



(37) 



(38) 



where we have used equations (I20t and ( I22> to recast (Ji?^ in terms of (JB^- Note that, neglecting the factor w against kjc^ in 
equation ( I33> . and therefore in equations ( I34> and ( I35t . effectively reduces to neglecting the term proportional to oj in equation 
( I19> . Thus, for the modes of interest, the condition cj^ ^ fcjc? is a statement about force balance in the vertical direction, which 
is made explicit in equation ( I38> . In this way, we can see how important perturbations in the density are, in the presence of strong 
toroidal fields (see also Balbus & Hawley 1991). For ua0 ^ Cg, even small variation in the toroidal component of the field can 
have an important impact on the dynamics of the perturbations. For this reason, the assumption of an incompressible MHD flow 
is not valid, whenever superthermal toroidal fields are considered. Note that, in order to recover the incompressible MRl when 
Ei = 0, for 1 = 1,2, 3, 4, we have not neglected the factor unity against {va^/cs)^, in equation ( I37> . 
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Fig. 4. — Left panels: Solutions to the full dispersion relation 1251 . Central panels: Solutions to the 4' order, approximate dispersion relation 1411 . with 
£2 = £3 = ^4 = 1. Right panels: Solutions to the 2""^ order, approximate dispersion relation 1441 . with £2 = ^3 = 1 and e4 = 0. All solutions correspond 
to a Keplerian disk with Cs = 0.05, vp^^ = 0.01, and vj^^ = 0.40. Open circles in upper panels indicate unstable modes. Note that the phase velocities of the 
two instabilities (seen in either the leftmost or central upper panels and corresponding to Region II and III in i|4. U are similar to the phase velocities, positive and 
negative respectively, of the slow mode (point-dashed line) in the limit of no rotation. The fast magnetosonic modes can barely be seen close to the left axis in 
the upper left panel. 



We now have all the elements to write equations ( I17> and ( I18> in terms of 5Br and JB^. Using equations ( I36> -( I38> . valid in 
the Umit lu"^ <C fcrCg, we obtain, in terms of the dimensionless variables. 



u'^SBr + 2iuj 



-(^y 



5B. 



oj'^SBa 



-(^) 



lUJ 



£4-; 



dlnn 
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dlnr 
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(kzVAz 



264 



UJ 



-'--m 



WA0 

kz VAz 
2' 



6Br 



SB. 



^{kzVAz)'^SB^ + ikzVA^VAzesSBr 



(39) 



(40) 



These equations are the generalization of the set of equations used to illustrate the physics behind the weak-field MRI as a system 
of masses coupled by a spring in a differentially rotating background. Indeed, in the incompressible limit and neglecting the 
curvature terms proportional to e,;, fori = 1, 2, 3,4, we recover the set of equations presented elsewhere (Balbus & Hawley 1992, 
1998). 

Setting the determinant of the linear system ( I39> - (I40> equal to zero and taking the limit ua^ ^ Cg provides the following 
approximate dispersion relation that is valid for strong toroidal fields^. 



- (k^ -I- k^VA^ + e2€iVA^)uj'^ - 2kzVA4,VAz{e2 + £3)1^ + klvA^ 



^A0 



2„,2 



-^ kivi 



dlnr 



e2e3WA0 



= 



(41) 



Note that we have not neglected the factor c^/v\j^ in the last term in equation ( I41> because its contribution is non-negligible 
at large wavenumbers. The solutions to the dispersion relation ( I4U . for a Keplerian disk with Cg = 0.05, vaz = 0.01, and 



' Note that, had we taken the opposite limit, i.e., Cs ^ ^A(^> we would have recovered the dispersion relation 1281 in the limit kr/k^ 
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Fig. 5. — The black dots represent unstable modes satisfying the approximate instability criteria 1431 . described in i|4.2l The dashed lines, labeled by fcj^ 
and k'^, are the limits of Regions I and II obtained analytically, also in i|4.2l The onset of instability III is labeled by k'^^ . As in Figure|3] we have assumed a 
Keplerian disk with Cs = 0.05 and i^Az = 0.01. 



VAip = 0.4, are shown in the central panels in Figure|3 For the sake of comparison, the left panels in the same figure show the 
solutions of the full dispersion relation ( I25> . The solutions to the approximate dispersion relation ( I41> are in excellent agreement 
with the solutions to the general dispersion relation j25l l for which lo^ ^ ^Ic^. 

Note that the term proportional to ei is not present in equation ( I4U . This feature has important conseq uences for us to 
understand the physics behind the stability of strongly magnetized compressible flows. It has been suggested ("Cur rv & Pudritd 
[1995) that the magnetic tension term B^6B^/rQ (i.e., the one proportional to ei in eq. IITI ) is responsible for the stabilization of 
long-wavelength perturbations via the restoring forces provided by strong toroidal field lines in incompressible MHD flows. This 
argument sounds compelling, but we can see from the last term in equation ( I39l l that the term proportional to ei is not dynamically 
important for compressible flows in which wa^ ^ Cg. At least in the radial direction, it is rather the term proportional to €2 the 
one governing the deviation compared to the stability properties of weak toroidal fields. This is in complete agreement with 
equation ( 13 8> . 



The dispersion relation (I41> is of the form 



62^^ + hiio + 60 = 



(42) 



For this 4**^ order equation to have complex roots (corresponding to unstable modes), its discriminant has to be negative, i.e.. 



Di{vA^, k.VAz) = -46^&I - 276^ + I660&2 - I2862&0 + 14462&f&o + 2566;] < 



(43) 



The modes satisfying this condition are shown as black dots in Figure|5] This analytical criterion agrees well with the numerical 
results for most of the parameter space (vAdJcs, kzVAz/cs) with the exception of some of the unstable modes close to the 
separatrix of the Regions 1 and 11, defined in q4. II 

Limiting wavenumbers for Regions I and II. — The modes satisfying the condition D^ = correspond to the limits of Regions I, 
II, and III in Figure|5l Their analytical expressions, however, are complicated. Yet, some more progress can be made by realizing 
that the solutions to the second order equation obtained by simply dropping the lo'^ term in equation ( 14 1> . 



(k^ + klv\^ + e2tiv\^)uj'^ + 2kzVA4,VAzie2 + £3)^ - k^vj^^^ 
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, 2 2 
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dlnn 
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dlnr 



<^2e3VA4, 



= 0, 



(44) 



constitute a very good approximation to the solutions of the dispersion relation ( I25> whenever the frequencies of the modes satisfy 
w^ <C 1. This can be appreciated by comparing the left and right panels in Figure |3 The physics behind this approximation is 
not as direct as the physics behind the condition cj^ <C k^cl, but it can also be understood in terms of force balance, this time 
in the radial direction. The dispersion relation ( I44> can be obtained by neglecting the term proportional to uj'^ in equation ( I39l l, 
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setting to zero the determinant of the resulting linear system given by equations (I39>-( I40> and taking the limit wa0 ^ Cg. This 
approximation is equivalent to neglecting the term proportional to a; in equation ( I17> and hence related to neglecting the radial 
acceleration experienced by a displaced fluid element. 

Setting the discriminant of equation ( I44> to zero, gives an equation in k^ with solutions that are the limiting wavenumbers for 
the onset of instabilities I and II in Figure|5j i.e.. 



D2ivA4>,kzVAz) = {kzVAz)'^ + 



dlnQ 
I 

dlnr 



^'A0 



"A<f, 
^2 



£4- 



^A0^ 



(kzVAz)'^ 



dinn 

> 

dlnr 



^M 



•^A^ 



^4 7r 



A0 



e4WA0 



0. 



(45) 



Here, we have set £2 = £3 = 1 but have explicitly left 64 to show that its contribution to the onset of instabilities I and II is not 
important when VA4, S> Cg, as long as we are considering a rotationally supported disk. We mention, however, that the numerical 
solutions show that the contribution of the term proportional to £4 is small but not negligible for the unstable modes in region III. 
Neglecting the terms proportional to 64, the solutions to equation (145 > are simply 



r.2 



iKvAzf = 

One of these solutions coincides always with /cbh (eq- 1271 ). 

{kfvAz) 



d\n^ 
) 

d In r 



± 



1 



"A4, 

r.2 



dhlVl 



d\nr 



and the other one is 



{kfvAzf = ^ 



(46) 



(47) 



(48) 



The modes with wavenumbers in the range [mm(fc^^, /c^^),m,aa:(fc^^, /c^^)] are unstable. In Figure |5] the critical curves 
k^_(vA4,) and ^^•^(ua^) are shown, with the proper normalization, as dashed lines. The critical wavenumber fc^^ in equation 



will be positive only for toroidal Alfven speeds larger than 



v\^ = V^ . (49) 

This is the critical value of the Alfven speed beyond which the modes with longest wavelength in Region I (see Fig. |5j begin 
to be stable. For a Keplerian disk, the epicyclic frequency coincides with the orbital frequency and thus, in dimensionless units, 
K^ = 1. In this case, the critical Alfven speed for fc^^ to be positive corresponds to VA4, = 0.223. This is the reason for which the 
long-wavelength modes are already stable in the second plot in the left panel in Figure[2 where VA4, = 0.25. 

Incidentally, we find that the values of toroidal Alfven speeds for which the standard MRJ gives the appropriate range of 
unstable modes are not as restricted to VA4, <C Cg but rather to VA4, <C y/ncs- For vacI> }t ^/^Cs, the standard MRI is stabilized 
at low wavenumbers. We point out that, Papaloizou & Szuszkiewicz (1992) found, by means of a global stability analysis of a 
compressible flow, that for a slim disk threaded only by a vertical field, the flow is stable if the vertical Alfven speed exceeds, 
within a factor of order unity, the geometrical mean of the sound speed and the rotational speed. In dimensionless units, this 
stability criterion translates into vaz ^ \/cI- 



The limiting case in which fc^^ — fc^^, is reached for 

v'l^ = 72^ . (50) 

Note that, for Cg = 0.05, this corresponds to a value for the critical toroidal Alfven speed of VAcf, = 0.316. This situation is to be 
compared with the mode structure in the third plot in the left panel in Figure^ where VA4, = 0.32. 

Limiting wavenumbers for Region III. — In the previous section we presented some useful analytical approximations to describe 
the dependence of the critical values of the toroidal Alfven speeds and wavenumbers defining Regions I and II on the different 
quantities characterizing the MHD flow. We could not, however, find simple analytical expressions to describe satisfactorily the 
corresponding behavior of the critical values defining Region III. We will describe next how the different unstable regions in 
Figure|3]depend on the magnitude of the sound speed and the steepness of the rotation profile. 

5. DISCUSSION 

In this section, we address several issues related to the importance of the curvature terms in determining the stability criteria 
obeyed by the solutions to the dispersion relation (I25> . We comment on some controversies raised by previous investigations 
that have treated the standard MRI taking into account, in various ways, either compressibility or curvature of the background 
magnetic field. We also comment on the importance in the outcome of the instabilities played by the magnetic tension produced 
b y toroidal filed lines in th e limit of cold MHD flows. We highlight the similarities and differences of our findings with the results 
of iCurrv & Pudritz ( 1995), who also found the emergence of a new (but different) instability for strong toroidal fields in the case 
of an incompressible MHD flow. Finally, we address the potential implications of our findings for shearing box models in which 
magnetic tension terms, induced by the curvature of the background field, are not considered. 
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we have considered a Keplerian disk with Cs = 0.05, haz = 0.01, and v/^^ = 0.4. 



5.1. Importance of Curvature Ternis 

In section fj^l'we mentioned that the terms proportional to ti, for i ~ 1, 2, 3, 4, are usually neglected in local stability analyses 
due to their 1/ro dependence. Some of these terms, ho'wever, are also proportional to the magnitude of the toroidal field. In 
this paper, ■we found that, ■when strong toroidal fields are considered, these terms led to substantial modifications to the stability 
criteria of MHD modes known to be valid in the limit of weak fields. After solving the full problem, we are in a better position 
to understand why this is the case. 

To illustrate the point, consider the ratio of the term proportional to fc^ to the one proportional to £2 in equation \\1\ and the 
ratio of the term proportional to fc^ to the one proportional to £3 in equation (I18> . i.e.. 



-Ri 



£2 "^A^ ^P 
ikz VAz SvAr 



and 



R2 



£3 VA4. 5vAr 



VAz 



(51) 



In order to ensure that the contributions due to curvature are negligible in a local analysis regardless of the magnitude of the 
toroidal field component, we should be able to ensure that the conditions i?i ^ 1 and i?2 <C 1 hold in the limit of large fc^ for 
any value of wa^ < 1- While it is encouraging that both dimensionless ratios are proportional to l/fc^, they are also proportional 
to the ratio of perturbed quantities, which we do not know a priori. It is only after having found the eigenfrequencies uj{kz) by 
taking into account all the curvature terms that we can properly address this issue. 

We can calculate how the ratios Ri and R2 depend on the wavenumber kz as follows. The ratio Ri can be recast using 
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Fig. 7. — The black dots represent unstable modes obtained from solving numerically the dispersion relation 1251 as a function of the toroidal Alfven speed. 
As an example, we have assumed Cs = 0.05 and ?jaz = 0.01. In each plot, we consider different values of the rotational profile, q = — dlnf2/dlnr. Note that, 
the highest value of the local toroidal Alfven speed considered here corresponds to the local circular velocity. 



equations dSSl i and ( l39b as 
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In a similar way, we can rewrite the ratio R2 using equation (I40> as 
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(53) 



For the sake of simplicity, let us consider a given value for the toroidal Alfven speed, e.g., VA4, = 0.4. Figure |6l shows the 
dependence of the ratios i?i and R2 on wavenumber for the unstable. The eigenfrequencies uj(kz) were obtained by solving 
equation ( I25> with e^ = 1, for i = 1,2, 3, 4, considering a Keplerian disk with Cg = 0.05 and vaz = 0.01. The ratios Ri and R2 
for the unstable modes (with wa^ — 8cs) in Regions II and III in Figure |3l are clearly identified. The complete mode structure 
corresponding to this case can be seen in the left panels of Figure|3 

It is important to stress that neither the real nor the imaginary parts of either i?i or R2 are negligible compared to unity even 
for Alfven speeds of order a few times the sound speed. In fact, for the unstable modes, the ratio Ri is of order unity and the ratio 
i?2 is in some cases larger than one by one order of magnitude. Their functional form is significantly different than the assumed 



5.2. Magnetorotational Instabilities with Superthermal Fields 

In we demonstrated that, when the toroidal magnetic field in a differentially rotating MHD flow becomes superthermal, 
three distinct instabilities can be identified, which we denote by roman numerals I, II, and III in Figure |3l We summarize the 
qualitative characteristics of these instabilities below. 
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Fig. 8. — The black dots represent unstable modes obtained from solving the dispersion relation 1251 numerically, for a Keplerian disk with vx^ = 0.01. In 
each plot, different values of the local sound speed, Cs, are considered. Note that in this case, the axes are not normalized by the pailicular value of the local 
sound speed, but rather by our initial choice of dimensionless variables, see J2| 



In contrast to the weak-field MRI, all three instabilities correspond to compressible MHD modes. Moreover, while the tra- 
ditional MRI corresponds to perturbations with negligible displacements along the vertical direction, this is not true for any of 
the three instabilities with superthermal toroidal fields. Instead, vertical displacements are an important characteristic of these 
instabilities and they occur with negligible acceleration, under a force balance between thermal and magnetic pressure. Finally, as 
in the case of the MRI, there is no significant acceleration along the radial direction but rather a force balance between magnetic 
tension, magnetic pressure, and thermal pressure. 

In Figures0and|8]we study numerically the dependences of the three instabilities on the properties of the background flow. As 
also shown in the case of the weak-field MRI (Balbus & Hawley, 1991), instability I occurs only in differentially rotating flows, 
with radially decreasing angular velocity. However, instability I also requires the presence of a non-negligible thermal pressure. 
Either a radially increasing angular velocity or a superthermal toroidal field can suppress instability I and hence the traditional 
MRI. 

Instability II is ubiquitous, whenever the background toroidal field of the flow is significantly superthermal. Indeed, it occurs 
even for flat (see Fig.Q or very cold (see Fig.|8|i flows. In a sense opposite to instability I, the steepness of the rotational profile 
determines the minimum unstable wavenumber, whereas the magnitude of the sound speed determines the minimum toroidal 
field strength required for the instability to o ccur. This inst ability seems to correspond to a generalization of the axisymmetric 
toroidal buoyancy (ATB) modes identified in iKim & Ostriker. C2Q0Q) . where the case Cs = was studied. In a similar way to 
instability II, the ATB modes with Cg = become unstable for all wavenumbers exceeding a critical value (for vertical modes 
this value is just given by /sbh)- When a finite sound speed is considered, however, thermal effects play an important role at small 
scales by completely stabilizing all the modes with wavenumbers larger than fc^^ (eq. 1481 ). 

Finally, instability III depends strongly on the rotational profile but very weakly on the sound speed. For rotationally supported 
flows (i.e., for ua^ ^ 1), instability III occurs only for significantly steep rotational profiles, e.g., q ~ \d\nVt/d\nr\ > 1.0, for 
the parameters depicted in FigureQ 
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Fig. 9. — The growth rate evolution of the different instabiUties defined in 34. ll for increasing magnetic field strength parameterized in terms of the ratio c^/v^ 
for a fixed pitch angle i = tan^^(i)Az/^A(/>) = 25°. Left panels: Growth rates of the unstable solutions to t he fu ll dispersion relation 1251 . when all curvature 
terms are taken into account. Right panels: Growth rates of the unstable solutions to the dispersion relation 1301 . i.e., when all curvature terms are neglected. 
In both cases, we have considered a vanishing ratio kr/k^, v\^ = 0.05, and a Keplerian disk. The dotted lines, in both panels, show the stabilization of the 
standard MRI as the magnetic field becomes superthermal (instability I). The soUd and dashed lines on the left panel show the growth rates corresponding to 
instability III and II respectively. As discussed in §3] these instabilities do not have a counterpart when the magnetic tension induced by bending of toroidal field 
lines is neglected. For completeness we have included, in the right panel, with solid and dashed lines the unstable solutions to equation (44) in Kim & Ostriker 
(2000). Instability II corresponds to a generalization of the ATB mode when thermal effects are accounted for 



5.3. Comparison to Previous Analytical Studies 

Soon after the original paper by Balbus & Hawlev ( 1991), Knobloch ( 1992) critiqued their approach to the study of local 
instabilities for lacking the contributions of curvature terms. iKnobloch (11992 ) formulated the stability analysis of a vertically 
unstratified, incompressible disk as an eigenvalue problem in the radial coordinate. He found that the presence of a toroidal 
field component changes the conditions for the presence of the instability as well as the character of the unstable modes from 
purely exponentials to overstable (i.e., i?e[tj] ^ 0). Ga mmie & Balbus (1994) argued against Knobloch's findings regarding 
overstability, stating that it arose as a consequence of having kept only small order terms (like wa/cs and wa/^o''o)- They 
concluded that these contributions would have been negligible had the flow been considered compressible. 



As we comment in ti4.2l Knobloch's dispersion relation is correct even in the limit Cg ^ wa0 (i-^-j without the necessity of 
imposing strict incompressibility). Formally speaking, the linear term in lo in equation ( I28> does break the symmetry of the 
problem allowing for unstable modes with Re{Lo) ^ 0. But it is also the case that, in the limit Cg ^ wa0, because of the relative 
magnitude of the coefficients in the dispersion relation (I28> . we do n ot ex pect the stability properties of the flow to differ greatly 
from those described by the incompressible MRI. As we mention in il2.3l in order to see significant differences, the Alfven speed 
would have to be of the order of the circular speed and therefore we do not expect the curvature terms to play a significant role 
on the stability of incompressible, rotationally supported flows. On the other hand, if we allow the MHD fluid to be compressible 
and consider the curvature of the background flow, the mode structure can be radically different from what is expected for the 
compressible MRI (c.f. Blaes & Balbus 1994). This is the case, even if the toroidal Alfven speed exceeds the sound speed by a 
factor of a few without the necessity of violating the condition of a rotationally supported disk (see Fig. Q. 

The stability of axisymmetric perturbations in weakly ionized and weakly magnetized shear flows was considered by 
iBlaes & Balbus ( 1994). They showed that, when ionization equilibrium is considered in the two-fluid approach, strong toroidal 
fields can fully stabilize the flow. As part of their study, they relaxed the Boussinesq approximation in the case of a single 
flu id and argued that, to a ll orders in the field strength, the magnitude of B^ does not affect the stability criterion. As noted 
bv lCurrv & Pudrita (119951) . this conclusion was reached because the terms proportional to B^/ro were not included in the local 
analysis. 

The behavior of the MRI in cold MHD shearing flows, has been addressed by Kim & Ostriker (2000). When performing their 
local analysis, these authors obtained the compressible version of the standard dispersion relation for the MRI and studied its 
solutions for different values of the ratio c^/v]^. Analyzing their dispersion relation (i.e., their equation [57] which is equivalent 
to equation 1301 in this study), Kim & Ostriker (2000) concluded that, when the magnetic field is superthermal, the inclusion 
of a toroidal component suppresses the growth rate of the MRI. Moreover, they found that, for a Keplerian rotation law, no 
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axisymmetric MRI takes place in very cold MHD flows if i < 30°, where i is the pitch angle of the local magnetic fields, 

i = tan"^(wAz/wA0)- 
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Fig. 10. — The real parts of the solutions to the dispersion relation 1251 corresponding to c? /iij = 0.01 for a pitch angle j = tan ^{vAz/^Aip) = 25°, 
a vanishing ratio kr/kz, s^Az = 0.05, and a Keplerian d isk. Left panel: Solutions to the full dispersion relation I25i . when all curvature terms are taken into 
account. Right panel: Solutions to the dispersion relation 1261 , i.e., when all curvature terms are neglected. Open circles indicate unstable modes. In both cases, 
the unstable modes with the shortest wavelength correspond to fcen^'Az- 



In 21 we showed that, depending on the strength of the toroidal field component, accounting for the finite curvature of the 
background magnetic field and the finite compressibility of the flow could be crucial in establishing which modes are subject 
to instabilities. In particular, we stated that both effects should be considered simultaneously whenever the local value of the 
toroidal Alfven speed exceeds the geometric mean of the local sound speed and the local rotational speed (for a Keplerian disk). 
However, this analytic criterion was found to be relevant for the modes with frequencies satisfying the condition lj^ <C klc^. 
Therefore, it is not obvious that we can trust this criterion in the limit Cs -^ 0. 

In order to see whether finite curvature effects do play a role in the stability of cold MHD flows we solved the complete 
dispersion relation ( 125 > for a pitch angle i — 25° and considered different values for the ratio c^/v\. Figure|5]shows the growth 
rates for the unstable modes of our dispersion relation (i.e., eq. 1251 ) and compares them to the ones of the dispersion relation 
obtained when the curvature terms are neglected (i.e., eq. 1301 ). In both cases, we have considered k^ S> fcr, waz = 0.05, and a 
Keplerian disk. The stabilization of the standard MRI (i.e., instability I) as the magnetic field becomes superthermal (wa > Cs) 
is evident (dotted lines in both panels). When the effects of magnetic tension are considered, not only does the growth rate of 
the MRI decrease faster for low values of c^lv\ but the modes with longest wavelengths are no longer unstable (e.g. when 
c^/v\ = 0.01). Because of this, the MRI is completely stabilized even for finite values of c^/v\. In contrast, when the curvature 
of the field lines is neglected, the growth rates decrease but the range of unstable modes remains unchanged as Cg ^ (right 
panel in Figure|9j; it is only when Cg = that the MRI is completely suppressed. For completeness, we present in FigureFTOIthe 
real parts of the solutions to the dispersion relation (125 > for i = 25°, kz ^ k^, vaz = 0.05, c^lv\ = 0.01 and a Keplerian disk, 
in the cases where curvature terms are considered (left panel) and neglected (right panel). The stabilization of the MRI at low 
wavenumbers and the emergence of instability III are evident. Note that, from Figure^] it is clear that the inclusion of magnetic 
tension terms can cause modifications to the mode structure when va ^ Cg, even when the toroidal and vertical components of 
the magnetic field are comparable. 



In g5j2| we mentioned that instability II seems to be a generalization of the axisymmetric toroidal buoyancy (ATB) modes, 
identified by Kim & Ostriker (2000), that accounts for finite temperature effects. Further indication that this is indeed the 
case can be found in Figure |9] where we have plotted (dashed lines) the growth rates corresponding to instability II and the 
one corresponding to the ATB mode (solutions of equation [44] with w^ <C k1v\^ in Kim & Ostriker 2000). Although finite 
compressibility suppresses instability II at large wavenumbers, it is clear that, as Cg -^ 0, the growth rates associated with 
instability II tend continuously to the growth rate of the cold ATB mode. For completeness, we have also included, in the right 
panel of Figure |5] the growth rate corresponding to the remaining unstable solution of equation (44) in Kim & Ostriker (2000) 
(solid line). This growth rate should be interpreted with great care since the aforementioned equation was derived under the 
condition lo^ ^ k1v\, which is not satisfied by the corresponding unstable mode. Nonetheless, we have included it to show the 
similarities that it shares with the growth rates corresponding to instability III (solid lines in the left panel) as Cg ^ . Note 
that the growth rates corresponding to instability III increase as Cg ^ and they saturate at Cg — 0. Although the higher critical 
wavenumber and the growth rate around this critical wavenumber seem to be the same for both instabilities, the differences 
between them at low wavenumbers is also evident. These differences become more dramatic as the pitch angle increases. 
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Curvature terms cannot, of course, be neglected in global treatments of magnetized accretion disks. It is, therefore, not 
surprising that new insta bilities, distinct from th e MRI, have akeady been found in global studies in which strong fields were 
considered. In particular. lCurrv & PudritzUlQQSl) performed a global stability analysis to linear axisymmetric perturbations of an 
incompressible, differentially rotating fluid, threaded by vertical and toroidal fields. They considered power-law radial profiles 
for the angular velocity and the toroidal and vertical components of the field. Each of these were parameterized as fi ex r^", 
B^ oc r~''+^, and B^ oc 7-^'^+^, respectively. Most of their analysis dealt with a constant vertical field and they allowed 
variations of the exponents (a, b), with the restriction that they correspond to a physical equilibrium state with a stationary 
pressure distribution. Although the majority of that paper dealt with global characteristics, they also performed a WKB analysis 
and concluded that, for 3/2 < a = b < 2 and ua0 < 1, the growth rate of unstable modes is suppressed on both short and long 
wavelengths and it approaches zero when VA<f, — > 1. On the other hand, for a = b y^ 2 and ua^ > 1, they found a new instability, 
with a growth rate that increased with wa0. They call this the Large Field Instability (LFI) and showed that it can be stabilized 
for sufficiently large vaz ■ 

It is worth mentioning the major qualitative differences between the LFI and the new instability discussed in fj3]that emerges 
for kz < /cbh after the stabilization of the MRI. Although it is true that, for our instability to be present, it is necessary for the 
toroidal Alfven speed to exceed the local sound speed, there is no need to invoke Alfven speeds larger than the local rotational 
speed. This is in sharp contrast with the LFI which only appears for va^, > 1. Regarding the range of unstable wavenumbers, the 
LFI remains unstable for k^ — > 0, albeit with diminishing growth rate for large values of vaz ■ This is not the case for the new 
instability present at low wavenumbers in our study. Th is can b e seen, for example, in the left panels in Figure 0] Perhaps the 
most noticeable difference is that the two instabilities in lCurrv & Pudritz ( 1995) that are present in the case a = 6 7^ 2 do not 
seem to coexist under any particular circumstances. The instability present for vA(j> < 1 reaches zero growth for wa^ -^ 1, while 
the LFI appears for vacj, = 1 and its growth rate is proportional to VA,p- When compressibility is considered, however, the two 
new instabilities found in our study can coexist even for Alfven speeds smaller than the local rotational speed. 

5.4. Implications for Shearing Box Simulations 

In an attempt to capture the most relevant physics without all the complexities involved in global simulations, the shearing 
box approach has been widely used in numerical studies of magnetized accretion disks (see, e.g., Hawley, Gammie, & Balbus 
1994). The aim of the shearing box approximation is to mimic a small region of a larger disk. The size of the box is usually 
Hz X 2ttHz X Hz, with Hz the thermal scale height of the isothermal disk. In this approach, it is common to adopt a pseudo- 
Cartesian local system centered at tq and in corotation with the disk with an angular frequency Qq, with coordinates x ^ r — tq, 
y — ro{(j) — flat), and z. The effects of differential rotation are then considered by imposing a velocity gradient in the radial 
direction. For a Keplerian accretion disk this is achieved by setting Vy = —{3/2)ilox. 

In most studies of unstratified shearing boxes, Alfven speeds rarely exceed the value of the local sound speed (see, e.g., Hawley, 
Gammie, & Balbus 19 95, 1 996). This is mainly because they are designed to simulate the mid-plane of the disk where the flow 
is relatively dense. In ii4.2l we have seen that, as long as the toroidal Alfven speed does not exceed the critical value u^ , ~ CgK, 
neglecting magnetic tension due to the curvature of toroidal field lines does not seem to affect the outcome of the MRI and hence 
the shearing box approach is well justified. However, when stratification is taken into account, usually by adopting a density 
profile of the form p ex exp[— z'^/(2iJ^)] in the case of isothermal disks, the steep drop in the density beyond a few scale heights 
can potentially lead to a magnetically dominated flow, with Alfven speeds larger than the critical value v^ ; = KCg. 

As discussed in the introduction. Miller & Stone (2000) carried out three-dimensional MHD simulations to study the evolution 
of a vertically stratified, isothermal, compressible, magnetized shear flow. The simulations were local in the plane of the disk 
but vertically extended up to ±5 thermal scale heights. This allowed them to follow the highly coupled dynamics of the weakly 
magnetized disk core and the rarefied magnetically-dominated (i.e., /3 < 1) corona that formed above the disk, for several (10 to 
50) orbital periods. 

Miller & Stone (2000) considered a variety of models, all sharing the same initial physical background, but different initial field 
configurations. In particular, they considered the following values: ilo = 10^^, 2Cg — 10"^ (so that Hz — ^/2cs/^o — 1), and 
tq = 100. We mention some of the results they obtained for the models with initial toroidal fields (BY), which were qualitatively 
similar to the zero net z-field (ZNZ) models. After a few orbital periods, the presence of a highly magnetized (with plasma 
/3 ~ 0.1 — 0.01) and rarefied (with densities two orders of magnitude lower than the disk mid-plane density) corona above ^ 2 
scale heights is evident. Within both the disk and the corona, the "toroidal" component of the field {By), favored by differential 
rotation, dominates the poloidal component of the field by more than one order of magnitude (with B^ ~ B^). 

We can compare the predictions of our study to the mode structure that one might expect from the standard compressible 
MRI for the particular values of sound and Alfven speeds found in the strongly magnetized corona by Miller & Stone (2000). 
To this end, we consider as typical (dimensionless) values Cg = 0.007, VA4, = 0.1 and vaz = 0.01, where we have assumed 
(3 = 2(cs/wa)^ — 0.01 and B^ = IQBz- The largest features that their simulations are able to accommodate are those with 
kz ~ 60 (corresponding to a wavelength of 10 in the vertical direction). As it can be seen in Figure^J the role of the curvature 
terms is not negligible in two different respects. First, it completely stabilizes the perturbations on the longest scales well inside 
the numerical domain. Second, the growth rate for the unstable modes is significantly reduced. 

It is difficult to extrapolate from the present work to addre ss how the instabilities discussed here would couple to buoyancy in 
the presence of a stratified medium like the one considered bv lMiller & Stone|(l2000h . Shearing boxes might also suffer from other 



22 



PES S AH & PSALTIS 






' / ' ' I : ' 
J ! 

Lji r 







1 



-2 



:^^^ 







curvature 



\ 



L 



I A.... ^ 

\ i --..... 

■ V ... I ........ . .*?».- 



.i 



i ! 

/ 



no curvature 



.1 



■ i I 



V . . . 





1 , , , , 1 , , , 


1 , , , , 


0.2 


curvature - 


0.1 


- 


0.0 


■ ^''' — ~^ 


i v„, ^^^y 


0.1 
0.2 







i 
i 


1 , , , , 


- 


no curvature 




i 




■ 


^^^^^^^ 


~~~~^^ 


"^^^-^'^ 


___^ 


- 


i 




i 
..... i ..... . 




- 



50 



100 

^zVaz 



150 



200 



50 



100 

^zVaz 



150 



200 



Fig. 1 1 . — The implication of our study for shearing box simulations. Lef t panels: Solutions to the full dispersion relation 1251 . when all curvature terms 
are taken into account. Right panels: Solutions to the dispersion relation 1261 . i.e., when all curvature terms are neglected. In both cases, we have considered 
Cs = 0.007, ?)Az = 0.01, vp^^ = 0.1, and a Keplerian disk. Open circles in upper panels indicate unstable modes. The vertical line indicates the minimum 
wavenumber (i.e., largest wavelength) that can be accommodated in the simulations of a strongly-magnetized corona above a weakly magnetized disk by 
iMiller & Stone. klOOO,) . 



problems when used to model strongly magnetized plasmas (e.g., the shearing sheet boundary conditions in the radial direction 
might not be appropriate for strong fields). The question is raised, however, about whether, because of their own Cartesian nature, 
they constitute a good approach at all to simulating compressible flows in which superthermal toroidal fields are present. Despite 
the fact that the generation of strongly magnetized regions via the MRI in stratified disks seems hard to avoid, their stability 
properties will ultimately depend on both the use of proper boundary conditions and proper accounting of the field geometry. 



6. SUMMARY AND CONCLUSIONS 

In this paper we have addressed the role of toroidal fields on the stability of local axisymmetric perturbations in compressible, 
differentially rotating, MHD flows, when the geometrical curvature of the background is taken into account. In order to accom- 
plish this task, without imposing restrictions on the strength of the background equilibrium field, we relaxed the Boussinesq 
approximation. In particular, we have studied under which circumstances the curvature terms, intimately linked to magnetic ten- 
sion in cylindrical coordinate systems, cannot be neglected. We have shown that the MRI is stabilized and two distinct instabilities 
appear for strong toroidal fields. At least for large wavenumbers, the structure of the modes seems to be the result of a purely 
local effect that is accounted for when compressibility and curvature terms are consistently taken into account. In particular, we 
have demonstrated that, even for rotationally supported cylindrical flows, both curvature terms and flow compressibility have to 
be considered if, locally, the toroidal Alfven speed exceeds the critical value given by u^ , ~ {K/n)cs^r (in physical units). 

There is little doubt that a realistic treatment of normal modes in magnetized accretion disks has to contemplate gradients in the 
flow variables over large scales and should, therefore, be global in nature. The results presented in this paper, however, provide 
the complete dispersion relation and, more importantly, analytic expressions for some of its solutions that should be recovered, 
in the appropriate limit, by a study of global modes in magnetized accretion disks, where compressibility effects are likely to be 
non negligible. We will address these issues in a future paper. 
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Fig. 12. — Solutions to the dispersion relation <251 for different values of the ratio k^/k^. The left, central, and right panels show the results for k^/kz = 
0.01, 0.1, and 1 respectively. In all three cases, we have considered Cg = 0.05, i;az = 0.01, vj^^ = 0.4, and a Keplerian disk. Open circles in upper panels 
indicate unstable modes. The same instabilities that were present in the case with vanishing ratio k^/k^ (left panels in Figure|4} can be clearly identified. 



APPENDIX A 

The central dispersion relation obtained in this study, i.e., equation ( I25> . was derived considering axisymmetric perturbations 
in both the vertical and the radial directions. Throughout the majority of our analysis, however, we focused our attention on the 
modes with vanishingly small ratios kr/k^. Particular emphasis has been given to the study of these modes in the literature of 
weakly mag netized, d ifferentially rotating disks, since these are the modes that exhibit the fastest growth rates (Balbus & Hawle^ 
fl992i.il998tlBalbusl2003.) . After having analyzed the role played by magnetic tension forces due to the finite curvature of strong 
toroidal field lines on the stability of these modes, we are in a better position to understand their effects on the modes for which 
the ratio kr/k^ is finite. 

In Figure [21 we present the solutions to the dispersion relation (125 > for three different ratios of the radial to the vertical 
wavenumber, i.e., k^/k^ — 0.01, 0.1, and 1. For the sake of comparison, we have used in this figure the same parameters that we 
used in obtaining the left panels in Figure^] for which the ratio k^/k^ was considered to be vanishingly small; we have assumed 
Cs — 0.05, vaz = 0.01, va0 = 0.4 and a Keplerian disk. In all the cases, the same instabilities (II and III) that were present in 
the case k^ ^ kr can be clearly identified. Note however, that even for very small ratios k^/k^ (e.g., left panels in Figure [T2t 
some of the modes that were stable in the case k^ ^ kr become unstable, albeit with negligible growth rate, when compared 
with the other unstable modes. It is also evident that in the limit k^ ^ fc^, the mode structure in Figure [21 tends continuously 
toward the mode structure in the left panels in Figurej^ It is this continuous behavior that ultimately justifies the study of modes 
with negligible ratio k^/k^ in a local stability analysis. 

In the case k,. — k^ (right panels in Figure [T2li . the value of the critical vertical wavenumbers for the onset of instabilities, 
i.e., fc^^ and k'^^ in equations ( I47> and ( 14 8 1 respectively, are different with respect to the case in which k^ S> fc^ by a factor 
\/2. This indicates that kz and kr play similar roles in establishing these critical wavenumbers. The growth rates of all these 
modes are reduced with respect to the case with k^ ^ kr- This behavior is similar to the one observed in the case of weak 
magnetic fields. It is important to stress that, even for the modes with comparable values of v ertical and radial wavenumbers, the 
general characteristics of the instabilities for strong toroidal fields that we discussed in !]4.1l are insensitive to the inclusion of a 
non-negligible kr- 

The completely new feature in Figure El is the appearance of another instability with a growth rate that does not seem to 
depended on wavenumber; the terms proportional to ikr in equation M5\ are crucial for the appearance of this new instability. 
The mode that is unstable seems to correspond to the mode that becomes the Alfven mode in the limit of no rotation. With 
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increasing values of the ratio kr/kz, the growth rates of the instabilities studied in J M.llgo to zero, but the new instability in 
Figure [T2lpersists . Note that for kr ~ fcz the growth rates of all the instabilities in Figure[21are comparable. 

As an aside, we point out that all the terms that are proportional to ikr, as opposed to fc^, in equation M5\ are also proportional 
to some factor e,; with i = 1, 2, 3, 4. All of these terms are negligible for sufficiently small ratios kr/k^ no matter how strong the 
toroidal field is. Indeed, if we consider perturbations with small enough radial wavelengths, at some point, curvature effects will 
not be important. However, this is not true in the vertical direction. In that case, we can ignore the curvature terms only when the 
toroidal field is weak. 



APPENDIX B 



In the various sections of the present study, we have seen that the toroidal component of the magnetic field does play a role 
in determining the stability criteria. In fact, for superthermal fields and quasi toroidal configurations, it dictates the values of 
some of the critical wavenumbers for the onset of instabilities (see ii4.2> . A particular, simple case, in which the importance of 
considering both compressibility and magnetic tension terms can be appreciated, is the study of modes with negligible frequency, 
i.e., with tj <C 1. We can obtain the wavenumber of these modes by imposing uj — to he a solution of the dispersion relation 
( I25> . We obtain, in physical dimensions. 
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In what follows, let us consider a rotationally supported disk whose rotational profile is not too steep (i.e., |d In il/d In 7-| is of 
order unity). The Alfven speed vaz appears in all three terms on the right hand side of equation (I54> and therefore it does not 
play a role in determining the relative magnitudes between them. Unlike the second term, the third term is no t necessarily small 
with respect to the first one, when superthermal fields are considered. In this case, it seems again (see M.2\ safe to neglect the 
curvature term proportional to ei6vA4, in equation ( I17> . However, had we neglected the curvature term proportional to e2Sp in 
equation (I17> or the one proportional to e^SvAr in equation ( I18> . we would have missed the important impact that the third term 
in equation (|53} has on the stability of modes with oj -^ 0, in the limit of strong toroidal fields (see Fig. 1^. This is somewhat 
counterintuitive because there does not seem to be any a priori indication about which of the magnetic tension terms (related to 
the curvature of the toroidal field component) is less relevant in the original set of equations (I16> -( I22> for the perturbations. This 
particular example illustrates the risks associated with neglecting terms that are not strictly 2"'* order in the perturbed quantities 
but rather address the geometrical characteristics of the background in which the (local) analysis is being carried out. 

From equation ( I54> it is also straightforward to see under which circumstances it is safe to neglect the curvature terms. For 
subthermal fields, fc° will not differ significantly from /cbh (eq. 1271 ) regardless of the geometry of the field configuration. This 
is because, if the field is weak enough {va <C Cg), no matter how strong of a (subthermal) i?^ component we consider, the 
second and third term are negligible with respect to the first one. We conclude this short analysis by commenting that, while 
a strong vertical field plays a stabilizing role, in the sense that it drives fc" toward small values leaving all modes with shorter 
wavelengths stable, the consequences of considering strong toroidal fields is a little more subtle as it can be seen in the evolution 
of the structure of the modes in Figure^ 
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